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Abstract. P-stars are compact stars made of up and down quarks in y6-equilibrium with 
electrons in a chromomagnetic condensate. We discuss p-stars endowed with super strong 
dipolar magnetic field which, following consolidated tradition in literature, are referred to 
as magnetars. We show that soft gamma-ray repeaters and anomalous X-ray pulsars can be 
understood within our theory. We find a well defined criterion to distinguish rotation pow- 
ered pulsars from magnetic powered pulsars. We show that glitches, that in our magnetars 
are triggered by magnetic dissipative effects in the inner core, explain both the quiescent 
emission and bursts in soft gamma-ray repeaters and anomalous X-ray pulsars. We account 
for the braking gUtch from SGR 1900-1-14 and the normal glitch from AXP IE 2259-1-586 
following a giant burst. We discuss and explain the observed anti correlation between hard- 
ness ratio and intensity. Within our magnetar theory we are able to account quantitatively 
for light curves for both gamma-ray repeaters and anomalous X-ray pulsars. In particular we 
explain the puzzling light curve after the June 18, 2002 giant burst from AXP IE 2259+586. 



Key words, pulsars: general - pulsars: individual: SGR 1900+14 - pulsars: individual: 
AXP IE 2259+586 

1. Introduction 

In few years since their discovery JHewish et al. 196St . pulsars have been identified with rotating 
neutron stars, first predicted theoretically by Baade & Zwicky (1934a, 1934b, 1934c), endowed 
with a strong magnetic field dPacini 1968IIGold 1968> . Nowadays, no one doubts that pulsars are 
neutron stars, even though it should be remembered that there may be other alternative explana- 
tions for pulsars. Up to present time it seems that there are no alternative models able to provide 
Send offprint requests to: Paolo Cea 
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as satisfactory an explanation for the wide variety of pulsar phenomena as those built around 
the rotating neutron star model. However, quite recently we have proposed (ICea 2004al b) a new 
class of compact stars, named p-stars, made of up and down quarks in y6-equilibrium with elec- 
trons in an abelian chromomagnetic condensate which is challenging the standard model based 
on neutron stars. 

In the present paper we investigate the properties of anomalous X-ray pulsars (AXPs) and soft 
gamma-ray repeaters (SGRs). For a recent review on the observational properties of anomalous 
X-ray pulsars see Mereghetti (1999), Mereghetti et al. (2002), Kaspi & Gavriil (2004), for soft 
gamma-ray repeaters see Hurley (1999), Woods (2003). Recently, these two groups have been 
linked by the discovery of persistent emission from soft gamma ray repeaters that is very similar 
to anomalous X-ray pulsars, and bursting activity in anomalous X-ray pulsars quite similar to 
soft gamma ray repeaters (see, for instance [Kaspi 2064| [Woods & Thompson 2004| l. Duncan & 
Thompson (1992) and Paczytiski (1992) have proposed that soft gamma-ray repeaters are pul- 
sars whose surface magnetic fields exceed the critical magnetic field Bqed - 4.4 lO'"* Gauss. 
Indeed, Duncan & Thompson (1992) refer to these pulsar as magnetars. In particular Duncan & 
Thompson (1995, 1996) argued that the soft gamma-ray repeater bursts and quiescent emission 
were powered by the decay of an ultra-high magnetic field. This interpretation is based on the 
observations that showed that these peculiar pulsars are slowing down rapidly, with an inferred 
magnetic dipole field much greater than the quantum critical field Bqeo, while producing steady 
emission at a rate far in excess of the rotational kinetic energy loss. 

The main purpose of this paper is to discuss in details p-stars endowed with super strong dipo- 
lar magnetic field which, following well consolidated tradition in literature, will be referred to 
as magnetars. We will show that soft gamma-ray repeaters and anomalous X-ray pulsars can 
be understood within our theory. The plan of the paper is as follows. In Section |5] we discuss 
the phenomenological evidence for the dependence of pulsar magnetic fields on the rotational 
period. We argue that there is a well defined criterion which allows us to distingue between ro- 
tation powered pulsars and magnetic powered pulsars. We explicitly explain why the recently 
discovered high magnetic field radio pulsars are indeed rotation powered pulsars. In Section |3 
we introduce the radio death Une, which in the P - P plane separated radio pulsars from radio 
quiet magnetic powered pulsars, and compare with available observational data. Sectionl^is de- 
voted to the glitch mechanism in our magnetars and their observational signatures. In Section RTTI 
we compare glitches in SGR 1900+14 and IE 2259+586, our prototypes for soft gamma-ray 
repeater and anomalous X-ray pulsar respectively. Sections f4.2l and l431 are devoted to explain the 
origin of the quiescent luminosity, the bursts activity and the emission spectrum during bursts. In 
Section lT^ we discuss the anti correlation between hardness ratio and intensity. In Section|5lwe 
develop a general formalism to cope with light curves for both giant and intermediate bursts. In 
Sections [5 . II through l5 .41 we careful compare our theory with the available light curves in litera- 
ture. In particular, we are able to account for the peculiar light curve following the June 18, 2002 
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giant burst from the anomalous X-ray pulsar IE 2259+586. Finally, we draw our conclusions in 
Section|6l 



2. ROTATION VERSUS MAGNETIC POWERED PULSARS 

As discussed in Cea (2004a) the structure of p-stars is determined once the equation of state ap- 
propriate for the description of deconfined quarks and gluons in a chromomagnetic condensate is 
specified. In general, the quark chemical potentials are smaller that the strength of the chromo- 
magnetic condensate. So that, up and down quarks occupy the lowest Landau levels. However, 
for certain values of the central energy density it happens that a fraction of down quarks must 
jump into higher Landau levels in the stellar core, leading to a central core with energy density 
Sc somewhat greater than the energy density outside the core. Now, these quarks in the inner core 
produce a vector cuiTent in response to the chromomagnetic condensate. Obviously, the quark 
current tends to screen the chromomagnetic condensate by a very tiny amount. However, since 
the down quark has an electric charge = -| e (e is the electric charge), the quark current 
generates in the core a uniform magnetic field parallel to the chromomagnetic condensate with 
strength Be - gf;; (here and in the following we shall adopt natural units /i - c - ks = 1). 
The inner core is characterized by huge conductivity, while outer core quarks are freezed into the 
lowest Landau levels. So that, due to the energy gap between the lowest Landau levels and the 
higher ones, the quarks outside the core cannot support any electrical current. As a consequence, 
the magnetic field in the region outside the core is dipolar leading to the surface magnetic field: 

- (If ' ^ 9^ ' 
R and Rc being the stellar and inner core radii respectively. In general the formation of the inner 
core denser than the outer core is contrasted by the centrifugal force produced by stellar rotation. 
Since the centrifugal force is proportional to the square of the stellar rotation frequency, this 
leads us to argue that the surface magnetic field strength is proportional to the square of the 
stellar period; 

Bs - Bi (t^) , (2) 
\ 1 sec/ 

where Bi is the surface magnetic field for pulsars with nominal period P = 1 sec. Indeed, 
in Fig. Q] we we display the surface magnetic field strength Bs inferred from (for instance, see 
[Manchester & Taylor 19771 1: 

Bs ^ 3.1 10''' 4FP Gauss , (3) 

versus the period. Remarkably, assuming Bi ^ 1.3 10'^ Gauss, we find the Eq. (|2j accounts 
rather well the infeiTed magnetic field for pulsars ranging from millisecond pulsars up to anoma- 
lous X-ray pulsars and soft-gamma repeaters. As a consequence of Eq. (|2j, we see that the dipolar 
magnetic field is time dependent. In fact, it is easy to find: 

Bsit) - Bo |1 + 2 ^ f ) , (4) 
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Fig. 1. Inferred magnetic field Bs plotted versus stellar period for 1 194 pulsars taken from the 
ATNF Pulsar Catalog jwww.atnf.csiro.au/research/pulsar/psrcal^ . Full line corresponds to Eq. (|2j 
withB, ^ 1.3 10^^ Gauss. 



where Bq indicates the magnetic field at the initial observation time. Note that Eq. Q implies 
that the magnetic field varies on a time scale given by the characteristic age. 
It is widely accepted that pulsar radio emission is powered by the rotational energy: 

1 



Er ^ - I a/ 
* 2 



so that, the spin-down power output is given by: 



- Er - - I cx) d) - A-n^ I — 

p3 



(5) 



(6) 



On the other hand, an important source of energy is provided by the magnetic field. Indeed, the 
classical energy stored into the magnetic field is: 



Assuming a dipolar magnetic field: 



B(r) = Bs 



or r > R 



Eq. leads to: 



Eb - - bI 
o 



(7) 



(8) 



(9) 



Now, from Eq. (0} the surface magnetic field is time dependent. So that, the magnetic power 
output is given by: 



Eb - - bI R^ - 



(10) 
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For rotation-powered pulsars it turns out that Eg <^ - Er. However, if the dipolar magnetic 
field is strong enough, then the magnetic power Eq. can be of the order, or even greater than 
the spin-down power. Thus, we may formulate a well defined criterion to distinguish rotation 
powered pulsars from magnetic powered pulsars. Indeed, until Eb < - Er there is enough 
rotation power to sustain the pulsar emission. On the other hand, when Eb > - Er all the 
rotation energy is stored into the increasing magnetic field and the pulsar emission is turned ofi^. 
In fact, in the next Section we will derive the radio death line, which is the line that in the P - P 
plane separates rotation-powered pulsars from magnetic-powered pulsars. In the remainder of 
this Section, we discuss the recently detected radio pulsars with very strong surface magnetic 
fields. We focus on the two radio pulsars with the strongest surface magnetic field: PSR J1718- 
3718 and PSR J1847-0130. These pulsars have inferred surface magnetic fields well above the 
quantum critical field Bqed above which some models ( [Baring & Harding \99%) predict that 
radio emission should not occur In particular, we have: 

P5/?71718-3718 dHobbs etal. 20041 P ^ 3.4 sec , Bs - 1 A 10^^ Gauss , 
f5/;/1847-0130 ( [McLaughlin et al. 2003> P - 6.7 sec , Bs ^ 9.4 10^^ Gauss . 

(11) 

Both pulsars have average radio luminosities and surface magnetic fields larger than that of AXP 
IE 2259-H586. Now, using Eqs. ([TO|l, (0 , together with / = § M/?^ and Eq. ((TT) we get: 

PSRjniS-371S - Er ^ 3.4 10^^ erg - , Eb ^ 4.7 10'^ erg - , 

^ ^ (12) 

PSRJIS47 -0130 - Er ^ 8.8 lO"^ erg ^ , Eb ^ 7.5 10^ erg ^ . 

We see that in any case: Eb < -Er, so that there is enough rotational energy to power the pulsar 
emission. 



3. RADIO DEATH LINE 

As discussed in previous Section, until Eb < - Er the rotation power loss sustains the pulsar 
emission. We have already shown that this explain the pulsar activity for pulsars with inferred 
magnetic fields well above the critical field Bqed- In this Section we explain why anomalous X- 
ray pulsars and soft gamma repeaters are radio quiet pulsars. When Eb ^ - Er all the rotation 
energy is stored into the increasing magnetic field and the pulsar emission is turned off. As a 
consequence pulsars with strong enough magnetic fields are radio quiet. Accordingly we see that 
the condition: 

Eb ^ - Er . (13) 

is able to distinguish rotation powered pulsars from magnetic powered pulsars. Now, using 
( [Manchester & Taylor 1977> 

i3IPP 

Bs ^ A > (14) 
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Fig. 2. Period derivative plotted versus stellar period for 1194 pulsars taken from the ATNF 
Pulsar Catalog ( |www.atnf.csiro.au/research/pulsar/psrcal| . Full line coiTesponds to Eq. ( I16t 



we recast Eq. ( fTsi into: 

P^P ^ \6n''R^ . (15) 
Using the canonical radius R ^ \OKm, we get: 

3 log(P) + log(P) ^ - 10 . (16) 

Equation ( I16> is a straight line, plotted in Fig. |2 in the log(P) - \og(P) plane. In 
Figure. |2 we have also displayed 1194 pulsars taken from the ATNF Pulsar Catalog 
www.atnf.csiro.au/research/pulsar/psrcat^ . We see that rotation powered pulsars, ranging from 
millisecond pulsars up to radio pulsars, do indeed he below our Eq. (I16> . Note that in Fig.|2]the 
recently detected high magnetic field pulsars are not included. However, we have already argued 
in previous Section that these pulsars have spin parameters which indicate that these pulsars are 
rotation powered. On the other hand, Fig.|2]shows that all soft gamma-ray repeaters and anoma- 
lous X-ray pulsars in the ATNF Pulsar Catalog he above our line Eq. M6\ . In particular, in Fig.|2l 
the pulsar above and nearest to the line Eq. (I16> corresponds to AXP IE 2259-1-586. So that, we 
see that our radio dead line, Eq. M6\ . correctly predicts that AXP IE 2259-1-586 is not a radio 
pulsar even though the magnetic field is lower than that in radio pulsars PSR J1718-3718 and 
PSR J1847-0130. We may conclude that pulsars above our dead line are magnetars, i.e. magnetic 
powered pulsars. The emission properties of magnetars are quite diff'erent from rotation powered 
pulsars: the emission from magnetars consists in thermal blackbody radiation form the surface. 
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4. GLITCHES IN MAGNETARS 

The origin of the dipolar magnetic field in p-stars is due to the inner core uniformly magnetized. 
In Figure |3] we display a schematic view of the interior of a p-star. The presence of the inner 




Fig. 3. Schematic view of the interior of a p-star Rc and R represent the inner core and stellar 
radii respectively; Sc is the energy density of the inner core, e is the energy density outside the 
core. 

core uniformly magnetized leads to well defined glitch mechanism in p-stars. Indeed, dissipative 
effects, which are more pronounced in young stars, tend to decrease the strength of the core mag- 
netic field. On the other hand, when decreases due to dissipation effects, then the magnetic 
flux locally decreases and, according to Lenz's law, induces a current which resists to the reduc- 
tion of the magnetic flux. This means that some quarks must flow into the core by jumping onto 
higher Landau levels. In other words, the core radius must increase. Moreover, due to very high 
conductivity of quarks in the core, we have: 

Be RI - constant , (17) 

which implies: 

6B,. SRc 

or 

^ . _ i ^ . (19) 
Rc 2 Be 

Equation J19> confirms that to the decrease of the core magnetic field, 6Bc < 0, it corresponds 
an increase of the inner core radius 6Rc > 0. This sudden variation of the radius of the inner core 
leads to glitches. Indeed, it is straightforward to show that the magnetic moment: 

m = Bs R^ = Be RI (20) 
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where we used Eq. ([TJ, must increase in the glitch. Using Eq. jl9L we get; 

— = ^ + 3-^--^>0, (21) 
m Be Rc Rc 

Another interesting consequence of the glitch is that the stellar radius R must decrease, i.e. the 
star contracts. This is an inevitable consequence of the increase of the inner core, which is char- 
acterized by an energy density higher then the outer core density e. As a consequence the 
variation of radius is negative: 6R < Q (see Fig. 13). In radio pulsar, where the magnetic energy 
can be neglected, conservation of the mass leads to: 



6R Ec-e /Ra^ 6Rc 

T ~ T~ \'r) 



(22) 



where e is the average energy density. In general, we may assume that is a constant of order 
unity. So that Eq. i22\ becomes: 

Note that the ratio can be estimate from Eq. Q once the surface magnetic field is known. 
We find that, even for magnetars, is of order 10"^ or less. So that our Eqs. Q and (I23> show 
that: 

6R 6Rc 

6R < , « — - . (24) 

R Rc 

As is well known, because the external magnetic braking torque, pulsars slow down according to 
(e.g. see Manchester & Taylor 1977): 

- V oc . (25) 

So that: 

6v 5m 61 6v 

— = 2 +3— . (26) 

V ml V 

From conservation of angular momentum we have: 

1^1 - \'-\ . (27) 

Moreover, from observational data it turns out that: 

6v 6v 

l-l » l-l , (28) 

V V 



6v ^ 6m ^ 5Rc 5 Be 

_ ~ 2 — =^2 — - ^ , (29) 

V m Rc Be 



so that Eq. ( I26> becomes: 



where we used Eqs. i21\ and (I19> . Equation i29\ does show that the variation of the radius of the 
inner core leads to a glitch. 

In rotation powered pulsar, starting from Eq. i22\ one can show that |^| ^ |^|. So that, taking 
into account Eq. J24> we recover the phenomenological relation Eq. ( I28t . A full account of 
gUtches in radio pulsar will be presented elsewhere. Glitches in magnetars are considered in the 
next Section, where we show that, indeed, Eqs. ( I24> and ( I28t hold even in magnetars. 
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The most dramatic effect induced by glitches in magnetars is the release of a huge amount of 
magnetic energy in the interior of the star and into the magnetosphere. To see this, let us consider 
the energy stored into the magnetic field in the interior of the magnetar We have: 

n2 



(30) 



where the first term is the energy stored into the core where the magnetic field is uniform. The 
variation of the magnetic energy Eq. i3Q\ caused by a glitch is easily evaluated. Taking into 
account Eq. iT% and <sc 1, we get: 

1 . (R? 6Rc 1 . . 6Rc 1 . J (R\^ 6Rc 

SE'"' ^ - - R^ \ — \ — R^ B^, — ^ - - R^ B^ \ — \ —. (31) 

^ 3 ^ \Rc Rc 2 Rc 3 ^ \Rc Rc 



Equation J31> shows that there is a decrease of the magnetic energy. So that after a glitch in 
magnetars a huge magnetic energy is released in the interior of the star We shall see that this 
energy is enough to sustain the quiescent emission. On the other hand, the glitch induces also a 
sudden variation of the magnetic energy stored into the magnetosphere. Indeed, from Eq. ([l) we 
find: 

^ = ^.3^-3^.^-3^.^>0. (32) 
Bs B(- Rc R Rc R Rc 

Thus, the magnetic energy stored into the magnetosphere: 

1 



^ -B's R' , (33) 



increases by: 



^ '-R^ BlC-^ .1'-^] ^ '-R^ Bl'-^ > . (34) 
3 ^ \ Bs 2 R j 3 ^ Bs 

This magnetic energy is directly injected into the magnetosphere, where it is dissipated by well 
defined physical mechanism discussed in Section 1431 and it is responsible for bursts in soft 
gamma-ray repeaters and anomalous X-ray pulsars. 

To summarize, in this Section we have found that dissipative phenomena in the inner core of a p- 
star tend to decrease the strength of the core magnetic field. This, in turn, results in an increase of 
the radius of the core 5Rc > 0, and in a contraction of the surface of the star, 5R < 0. We have also 
shown that the gUtch releases an amount of magnetic energy in the interior of the star and injects 
magnetic energy into the magnetosphere, where it is completely dissipated. Below we will argue 
that these magnetic glitches are responsible for the quiescent emission and bursts in gamma-ray 
repeaters and anomalous X-ray pulsars. Interestingly enough, in Cheng et al. (1995) it was shown 
that SGR events and earthquakes share several distinctive statistical properties, namely: power- 
law energy distributions, log-symmetric waiting time distributions, strong correlations between 
waiting times of successive events, and weak correlations between waiting times and intensities. 
These statistical properties of bursts can be easily understood if bursts originate by the release 
of a small amount of energy from a reservoir of stored energy. As a matter of fact, in our theory 
the burst activity is accounted for by the release of a tiny fraction of magnetic energy stored in 
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magnetars. Even for giant bursts we find that the released energy is a few per cent of the magnetic 
energy. Moreover, Hurley et al. (1994) noticed that there is a significantly statistical similarity 
between the bursts from SGR 1806-20 and microglitches observed from the Vela pulsar with 
1^1 ~ 10"''. So that we see that these early statistical studies of bursts are in agreement with 
our theory for bursts in magnetars. Even more, we shall show that after a giant glitch there is an 
intense burst activity quite similar to the setthng earthquakes following a strong earthquake. 

4.1. BRAKING GLITCHES 

In Sectionl^we found that magnetic glitches in p-stars lead to: 

- " - ^ > ■ (35) 
V Be 

Since there is variation of both the inner core and the stellar radius, the moment of inertia of 
the star undergoes a variation of 61. It is easy to see that the increase of the inner core leads 
to an increase of the moment of inertia /; on the other hand, the reduction of the stellar radius 
implies 61 < 0. In radio pulsar, where, by neglecting the variation of the magnetic energy, the 
conservation of the mass leads to Eq. ( I23t . one can show that: 

^ . ^ < . (36) 
I R 

Moreover, from conservation of angular momentum: 



61 _ _6v 
7 ~ ~ V 



(37) 



it follows: 



Sv 6R 6Be 6v 

0< - - - «-^ ^ _. (38) 

V R Br V 



For magnetars, namely p-stars with super strong magnetic field, the variation of magnetic energy 
cannot be longer neglected. In this case, since the magnetic energy decreases, we have that the 
surface contraction in magnetars is smaller than in radio pulsars. That means that Eq. J24> holds 
even for magnetars. Moreover, since in radio pulsars we known that: 

^ -^<10-, (39) 

V I R 

we see that in magnetars the following bound must hold: 

- f < 10-' . (40) 

As a consequence we may write: 

'I = I'l) , (E) (E) . £^ <o. 

/ \ / /core \ I I surf \ I I surf R 

As we show in a moment, the variation of the moment of inertia induced by the core is positive. 
So that if the core contribution overcomes the surface contribution to 61 we have a braking glitch 
where = ^ > 0. 

V P 



(41) 
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We believe that the most compelling evidence in support to our proposal comes from the anoma- 
lous X-ray pulsar AXP IE 2259+586. As reported in Woods et al. (2004), the timing data 
showed that a large glitch occurred in AXP IE 2259+586 coincident with the 2002 June gi- 
ant burst. Remarkably, at the time of the giant flare on 1998 August 27, the soft gamma ray 
repeater SGR 1900+14 displayed a discontinuous spin-down consistent with a braking glitch 
JWoods et al. 1999"^^ . Our theory is able to explain why AXP IE 2259+586 displayed a nor- 
mal glitch, while SGR 1900+14 sufi'ered a braking gUtch. To see this, we recall the spin-down 
parameters of these pulsars: 

SGR 1900 + 14 P ^ 5.16 sec, P ^ 1.1 10"'" , Bs ^ 7.4 lO"^ Gauss , 



AXP l£2259 + 586 P ^ 6.98 sec , P ^ 2.0 10"'^ , Bs ^ 1.2 lO'"* Gauss 



(42) 



For canonical magnetars with M ^ 1 .4 Mq and radius R ^ 10 Km, we have ■\JgH ^ 0.55 GeV. 
So that, using 1 GeV^ - 5.12 10^^ Gauss, we rewrite Eq. ([0 as: 

Bs ^ 1.54 lO"'^^) Gauss. (43) 

Combining Eqs. J42l i and J43> we get: 

(R \^ 

— 4.81 10"^ , 

^ 3 (44) 

AXP 1£2259 + 586 j ^ 0.78 10"^ . 

According to Eqs. J35> . (j^) and (I41> . to evaluate the sudden variation of the frequency and 
frequency derivative, we need 6R and dR,.- These parameters can be estimate from the energy 
released during the giant bursts. In the case of AXP IE 2259+586, the giant 2002 June burst fol- 
lowed an intense burst activity which lasted for about one year Woods et al. (2004), assuming a 
distance of 3 kpc to IE 2259+586, measured an energy release of 2.7 10^^ ergs and 2. 1 10"*' ergs 
for the fast and slow decay intervals, respectively. Due to this uncertainty, we conservatively 
estimate the energy released during the giant burst to be: 

AXP 1£2259 + 586 £w ^ 1.0 10"*° ergs . (45) 

On 1998 August 27, a giant burst from the soft gamma ray repeater SGR 1900+ 14 was recorded. 
The estimate energy released was: 

5GP 1900+ 14 Ebursr ^ 1.0 lO'*'* ergs . (46) 

As we have already discussed in Sect.|3 the energy released during a burst in magnetars is given 
by the magnetic energy directly injected and dissipated into the magnetosphere, Eq. i34l . We 
rewrite Eq. J34> as 

6ET . 1 p3 b2 ^ . 2.6 10^ ergs f ^ . (47) 

3 ^ Bs ^ UOi^GawiJ Bs 

So that, combining Eqs. j47l l. j46l l. j45l l and J42> we get: 



6Bs 6R,. , 

5GP 1900+ 14 ^ —L ^ 0.70 10^2 

Bs Rc 

6Bs SRc 9 

AXP 1£2259 + 586 ^ —L ^ Q21 10"^ 

Bs Rc 



(48) 
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Thus, according to Eq. ( I35> we may estimate the sudden variation of v: 

- - 2 — ^ ~ 10-2 , (49) 
V R,. 

for both glitches. On the other hand we have: 

(?) " T ^ (I)' ^ ' f (I)' ^ • <»> 

\ J /core S \ K / Kr L \ K ! Kr 



leading to: 



SGR19QQ+IA {^\ ^ 3.34 10-"^ 



AXP 1£2259 + 586 



(-1 

I / I CO 

(^) . 1.34 10-^ . 

\ I I core 



(51) 



On the other hand, we expect that during the giant glitch (y) ^ ~ 10"^. As a consequence, for 
AXP IE 2259+586 the core contribution is negligible with respect to the surface contribution to 
61. In other words, AXP IE 2259+586 displays a normal glitch with ^ ~ 10 On the contrary, 
Eq. (ED indicates that SGR 1900+14 suffered a braking glitch with ^ ~ -3.34 10 giving: 

AP ^ 1.72 10"^ sec . (52) 

We see that our theory is in agreement with observations, for a glitch of size ^ - 4.24(1 1) 10"^ 
was observed in AXP IE 2259+586 which preceded the burst activity d Woods et al. 2004t . 
Moreover, our theory predicts a sudden increase of the spin-down torque according to Eq. j49> . 
In Woods et al. (2004) it is pointed out that it was not possible to give a reliable estimate of 
the variation of the frequency derivative since the pulse profile was undergoing large changes, 
thus compromising the phase alignment with the pulse profile template. Indeed, as discussed in 
Sect. 15.11 soon after the giant burst AXP IE 2259+586 suffered an intense burst activity. Now, 
according to our theory, during the burst activity there is both a continuous injection of magnetic 
energy into the magnetosphere and variation of the magnetic torque explaining the anomalous 
timing noise observed in IE 2259+586. In addition. Woods et al. (1999b) reported a gradual 
increase of the nominal spin-down rate and a discontinuous spin down event associated with the 
1998 August 27 flare from SGR 1900+14. Extrapolating the long-term trends found before and 
after August 27, they found evidence of a braking glitch with AP ^ 0.57(2) 10"^ sec. In view 
of our theoretical uncertainties, the agreement with our Eq. (15 2> is rather good. To our knowl- 
edge the only attempt to explain within the standard model the braking glitch observed in SGR 
1900+14 is done in Thompson et al. (2000) where it is suggested that violent August 27 event 
involved a glitch. The magnitude of the ghtch was estimated by scahng to the largest glitches 
in young, active pulsars with similar spin-down ages and internal temperature. In this way they 
deduced the estimate |^| ~ 10"^ to 10 However, this explanation overlooks the well known 
fact that radio pulsars display normal glitches and no braking glitches. Moreover, there is no ar- 
guments to explain why AXP IE 2259+586 displayed a normal glitch instead of a braking glitch. 
Let us conclude this Section by briefly discussing the 2004 December 27 giant flare from SGR 
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1806-20. During this tremendous outburst SGR 1806-20 released a huge amount of energy, 
Eburst ~ 10^^ ergs. Using the spin-down parameters reported in Mereghetti et al. (2005a): 

SGR 1806-20 P ^ 7.55 sec, P ^ 5.5 10""' , Bs ^ 2.0 lO'^ Gauss , (53) 

we find: 

6Bs 6Rc , 
1806-20 ^ ^ — £ ^ 9.6 10 2 . (54) 

Bs Rc 

Thus, we predict that SGR 1806-20 should display a gigantic braking glitch with ^ ^ 2.4 10"^, 
or : 

^ 1.8 10"' sec . (55) 

4.2. QUIESCENT LUMINOSITY 

The basic mechanism to explain the quiescent X-ray emission in our magnetars is the internal 
dissipation of magnetic energy. Our mechanism is basically the same as in the standard magnetar 
model based on neutron star ( [Duncan & Thompson 1996) . Below we shall compare our proposal 
with the standard theory. In Section@]we showed that during a glitch there is a huge amount of 
magnetic energy released into the magnetar: 

As in previous Section, we use SGR 1900+14 and AXP IE 2259+586 as prototypes for soft 
gamma ray repeaters and anomalous X-ray pulsars, respectively. Using the results of Sect. 14.11 
we find: 

SGR 1900 + 14 - 5£™ ^ 3.0 lO'*^ erg — , 

(57) 

AXP 1£2259 + 586 - 6Ef ^ 4.8 10"^^ erg — . 

Rc 

This release of magnetic energy is dissipated leading to observable surface luminosity. To see 
this, we need a thermal evolution model which calculates the interior temperature distribution. 
In the case of neutron stars such a calculation has been performed in VanRiper (1991), where it 
is showed that the isothermal approximation is a rather good approximation in the range of inner 
temperatures of interest. The equation which determines the thermal history of a p-star has been 
discussed in Cea (2004a) in the isothermal approximation: 

dt 

where L,, is the neutrino luminosity, Ly is the photon luminosity and Cy is the specific heat. 
Assuming blackbody photon emission from the surface at an eff'ective surface temperature Ts 
we have: 

Ly = AjtR^o-sbT^s . (59) 
where os b is the S tefan - Boltzmann constant. In Cea (2004a) we assumed that the surface and 
interior temperature were related by: 

^ = 10"^ fl . (60) 



— = - (Lv + Ly) , (58) 
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Equation j60> is relevant for a p-star which is not bare, namely for p-stars which are endowed 
with a thin crust. The vacuum gap between the core and the crust, which is of order of several 
hundred /erm/s, leads to a strong suppression of the surface temperature with respect to the 
core temperature. The precise relation between Ts and T could be obtained by a careful study 
of the crust and core thermal interaction. In any case, our phenomenological relation Eq. ( I60t 
allows a wide variation of Ts , which encompasses the neutron star relation (see, for instance, 
IGundmundsson et al. 1983t . Moreover, our cooling curves display a rather weak dependence on 
the parameter a in Eq. i60\ . Since we are interested in the quiescent luminosity, we do not need 
to known the precise value of this parameter. So that, in the following we shall assume a ~ 1 . In 
other words, we assume: 

Ts ^ IQ-^ T . (61) 

Obviously, the parameter a is relevant to evaluate the surface temperature once the core temper- 
ature is given. Note that, in the relevant range of core temperature T ~ 10^ °K, our Eq. (16 1> is 
practically identical to the parametrization adopted in Duncan & Thompson (1996) within the 
standard magnetar model: 

T 



Ts ^ 1.3 10^ °K 



The neutrino luminosity Ly in Eq. (I58> is given by the direct yS-decay quark reactions, the domi- 
nant cooling processes by neutrino emission. From the results in Cea (2004a), we find: 

Ly ^ 1.22 10^^ — tI , (63) 
s 

where Tg is the temperature in units of 10'' °K. Note that the neutrino luminosity Ly has the 
same temperature dependence as the neutrino luminosity by modified URCA reactions in neutron 
stars, but it is more than two order of magnitude smaller. From the cooling curves reported in 
Cea (2004a) we infer that the surface and interior temperature are almost constant up to time r ~ 
10^ years. Observing that magnetars candidates are rather young pulsar with Tage ^ 10^ years, 
we may estimate the average surface luminosities as: 

Ly ^ ^ . (64) 

Tage 

We assume Tage ~ for SGR 1900+14. On the other hand, it is known that for AXP IE 
2259-1-586 Tage ~ 10^ years <K Tc. So that, we get: 

ii ere SR,- 

SGR 1900 +U Ly ^ 1.3 10^^ — — , 

(65) 

— /■' I S' iti\ ^ 

AXP IE2259 + 586 Ly ^ 1.5 10 



35 erg 6Rc 



s Rc 

We see that it is enough to assume that SGR 1900+14 suffered in the past a glitch with ^ ~ 10"^ 
to sustain the observed luminosity Ly ~ 10^^ ^ (assuming a distance of about 10 kpc). In the 
case of AXP IE 2259H-586, assuming a distance of about 3 kpc, the observed luminosity is 
Ly ~ 10^'^ — , SO that we infer that this pulsar had suffered in the past a giant gUtch with 
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~ 10 \ quite similar to the recent SGR 1806-20 glitch. 
Let us discuss the range of validity of our approximation. Equation ( I64l l is valid as long as Ly 
dominates over Ly, otherwise the star is efficiently cooled by neutrino emission and the surface 
luminosity saturates to L™"^. We may quite easily evaluate this limiting luminosity from L™°^ ^ 
Ly. Using Eq. J61t and R ^ iQ Km, we get: 

Ijnax ^ 4 2 10" ^ . (66) 

' s 

Note that, since our neutrino luminosity is reduced by more than two order of magnitude with 
respect to neutron stars, L™ * is about two order of magnitude greater than the maximum allowed 
surface luminosity in neutron stars ( |Van Riper 1991^ . Thus, our theory allows to account easily 
for observed luminosities up to 10^* 

Let us, finally, comment on the quiescent thermal spectrum in our theory. As already discussed, 
the origin of the quiescent emission is the huge release of magnetic energy in the interior of the 
magnetar Our previous estimate of the quiescent luminosities assumed that the interior tempera- 
ture distribution was uniform. However, due to the huge magnetic field, the thermal conductivity 
is enhanced along the magnetic field. This comes out to be the case for both electron and quarks, 
since we argued that the magnetic and chromomagnetic fields are aligned . As a consequence, we 
expect that the quiescent spectrum should be parameterized as two blackbodies with parameter 
R\ ,Ti and R2 , T2, respectively. Since the blackbody luminosities Lyi and Ly2 are naturally of the 
same order, our previous estimates for the quiescent luminosities are unaffected. Moreover, since 
the thermal conductivity is enhanced along the magnetic field, the high temperature blackbody, 
with temperature T2, originates from the heated polar magnetic cups. Thus we have: 

Ri < R , R2< I Km , 

(67) 

T2 > T\ , Rl > R2 , Lyl ^ Ly2 . 

Note that there is a natural anticorrelation between blackbody radii and temperatures. It is cus- 
tomary to fit the quiescent spectrum of anomalous X-ray pulsars and soft gamma ray repeaters 
in terms of blackbody plus power law. In particular, it is assumed that the power law compo- 
nent extends to energy greater than an arbitrary cutoff energy E cutoff - 2KeV. It is worthwhile 
to stress that these parameterizations of the quiescent spectra are in essence phenomenolog- 
ical fits. Indeed, within the standard magnetar model ( [Duncan & Thompson 1996) the power 
law should be related to hydromagnetic wind accelerated by Alfven waves. The luminosity of 
the wind emission should increase with magnetic field strength as L„,„£/ ^ Lpi ~ . On the 
other hand, the blackbody luminosities should scale as B^ '* ( [Duncan & Thompson 1996) . So that 
the ratio Lpi/LgB decreases with increasing magnetic field strengths, contrary to observations 
dMarsden & White 2001> . In our theory well defined physical arguments lead to the two black- 
body representation of the quiescent spectra, whose parameters are constrained by our Eq. (|67). 
As a matter of fact, we have checked in literature that the quiescent spectra of both anomalous 
X-ray pulsars and soft gamma ray repeaters could be accounted for by two balckbodies. For in- 
stance, in Morii et al. (2003) the quiescent spectrum of AXP IE 1841-045 is well fitted with the 
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standard power law plus blackbody {reduced /do f ^ 1.11), nevertheless the two blackbody 
model gives also a rather good fit (reduced ;t^^/c/o/ - 1.12). Interestingly enough the blackbody 
parameters: 

= 5.7 t?A Km , r, = 0.47 + 0.02 KeV , 

(68) 

R2 = 0.36!;j{J^ Km , T2 = 1.5 ^[i:^ KeV , 

are in agreement with Eq. ( I67t . Moreover, assuming that the power law component in the 
standard parametrization of quiescent spectra account for the hot blackbody component in our 
parametrization, we find that the suggestion Lyi ^ Ly2 in Eq. i61\ is in agreement with observa- 
tions JMarsden & White 200 1> . Another interesting consequence of the anisotropic distribution 
of the surface temperature due to strong magnetic fields is that the thermal surface blackbody 
radiation will be modulated by the stellar rotation. As a matter of fact, Ozel (2002) argued that 
the observed properties of anomalous X-ray pulsars can be accounted for by magnetars with a 
single hot region. It is remarkable that our interpretation explains naturally the observed change 
in pulse profile of SGR 1900-1-14 following the 1998 August 27 giant flare. It seem that our pic- 
ture is in fair qualitative agreement with several observations. However, any further discussion 
of this matter goes beyond the aim of the present paper. 

4.3. BURSTS 

In the present Section we discuss how glitches in our magnetars give rise to the burst activity 
from soft gamma-ray repeaters and anomalous X-ray pulsars. We said in Sect.|4lthat the energy 
released during a burst in a magnetar is given by the magnetic energy directly injected into 
the magnetosphere, Eq. (I34> . Before addressing the problem of the dissipation of this magnetic 
energy in the magnetosphere, let us discuss what are the observational signatures at the onset 
of the burst. Observations indicate that at the onset of giant bursts the flux displays a spike 
with a very short rise time fi followed by a rapid but more gradual decay time t2. According 
to our previous discussion, the onset of bursts is due to the positive variation of the surface 
magnetic field dBs , which in turn implies an sudden increase of the magnetic energy stored in 
the magnetosphere. Now, according to Eq. ( 0^ we see that almost all the magnetic energy is 
stored in the region: 

R < r < IQR . (69) 

So that the rise time is essentially the time needed to propagate in the magnetosphere the infor- 
mation that the surface magnetic field is varied. Then we are lead to: 

fi 9 ^ 3 10"^ sec , (70) 

which indeed is in agreement with observations. On the other hand, in our proposal the decay 
time f2 depends on the physical properties of the magnetosphere. It is natural to identify t2 with 
the time needed to the system to react to the sudden variation of the magnetic field. In other 
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words, we may consider the magnetosphere as a huge electric circuit which is subject to a sudden 
increase of power from some external device. The electric circuit reacts to the external injection 
of energy within a transient time. So that, in our case the time t2 is a function of the geometry and 
the conducting properties of the magnetosphere. In general, it is natural to expect that t\ t2 
and the time extension of the initial spike is: 

S tspike - h - h - h ■ (71) 

Remarkably, observations shows that the observed giant bursts are characterized by almost the 
same 6 tspike'- 

5 tspike - t2 ^ 0.1 sec , (72) 

signalling that the structure of the magnetosphere of soft gamma-ray repeaters and anomalous 
X-ray pulsars are very similar. Since the magnetic field is varied by 6Bs in a time dtspike, then 
from Maxwell equations it follows that it must be an induced electric field. To see this, let us 
consider the dipolar magnetic field in polar coordinate: 



(73) 



Bs R' sine 
Be T , 

0, 

Thus, observing that is the time derivative of the magnetic field it is easy to find the induced 
azimuthal electric field: 

SBs sine 

= +— ^ — , r > R . (74) 

^ispike ^ 

To discuss the physical eff^ects of the induced azimuthal electric field Eq. J74> . it is convenient 
to work in the co-rotating frame of the star We assume that the magnetosphere contains a neu- 
tral plasma. Thus, we see that charges are suddenly accelerated by the huge induced azimuthal 
electric field and thereby acquire an azimuthal velocity ^ 1 which is directed along the 
electric field for positive charges and in the opposite direction for negative charges. Now, it is 
well known that relativistic charged particles moving in the magnetic field B(r), Eq. il3\ . will 
emit synchrotron radiation (see for instance (Wallace 19771 [Ginzburg 1979| . As we discuss be- 
low, these processes are able to completely dissipate the whole magnetic energy injected into the 
magnetosphere following a glitch. However, before discuss this last point in details, we would 
like to point out some general consequences which lead to well defined observational features. 
As we said before, charges are accelerated by the electric field thereby acquiring a relativistic 
azimuthal velocity. As a consequence, they are subject to the drift Lorentz force F - q\^x B(r), 
whose radial component is: 

(R\^ 

F, - -qv^Bg ^ +qVy,Bs sine I— j , (75) 

while the 6 component is: 



(R \^ 
—j ■ (76) 
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The radial component pushes both positive and negative charges radially outward. Then, we 
see that the plasma in the outermost part of the magnetosphere is subject to a intense radial cen- 
trifugal force, so that the plasma must flow radially outward giving rise to a blast wave. On the 
other hand, Fg is centripetal in the upper hemisphere and centrifugal in the lower hemisphere. 
As a consequence, in the lower hemisphere charges are pushed towards the magnetic equato- 
rial plane cos = 0, while in the upper hemisphere (the north magnetic pole) the centripetal 
force gives rise to a rather narrow jet along the magnetic axis. As a consequence, at the onset 
of the giant burst there is an almost spherically symmetric outflow from the pulsar together with 
a collimated jet from the north magnetic pole. Interestingly enough, a fading radio source has 
been seen from SGR 1900-H14 following flie August 27 1998 giant flare (IFrail et al. 1999> . The 
radio afterglow is consistent with an outflow expanding subrelativistically into the surrounding 
medium. This is in agreement with our model once one takes into account that the azimuthal 
electric field is rapidly decreasing with the distance from the star, so that < 1 for the plasma in 
the outer region of the magnetosphere. However, we believe that the most compelling evidence in 
favour of our proposal comes from the detected radio afterglow following the 27 December 2004 
gigantic flare from SGR 1806-20 JGaensler et al. 20051 ICameron et al. 20051 |Wang et al. 2005| 
IGranot et al. 20051 iGelland et al. 20051 . Indeed, the fading radio source from SGR 1806-20 has 
similar properties as that observed from SGR 1900+14, but much higher energy. The interesting 
aspect is that in this case the spectra of the radio afterglow showed clearly the presence of the 
expected spherical non relativistic expansion together with a sideways expansion of a jetted ex- 
plosion (see Fig. 1 in lGaensler et al. 2005l and Fig. 1 in lCameron et al. 2005> . Finally, the lower 
limit of the outflow E > 10^'*^ ergs JGelfand et al. 20051 implies that the blast wave and the jet 
dissipate only a small fraction of the burst energy which is about lO"*^ ergs (see Section ITTl . 
Thus, we infer that almost all the burst energy must be dissipated into the magnetosphere. In 
the co-rotating frame of the star the plasma, at rest before the onset of the burst, is suddenly 
accelerated by the induced electric field thereby acquiring an azimuthal velocity ^ 1 . Now, 
relativistic charges are moving in the dipolar magnetic field of the pulsar. So that, they will lose 
energy by emitting synchrotron radiation until they come at rest. Of course, this process, which 
involves charges that are distributed in the whole magnetosphere, will last for a time f^,-, much 
longer that Stspike- Actually, tdis will depend on the injected energy, the plasma distribution and 
the magnetic field strength. Moreover, one should also take care of repeated charge and photon 
scatterings. So that it is not easy to estimate f^,-, without a precise knowledge of the pulsar magne- 
tosphere. At the same time, the fading of the luminosity with time, the so-called light curve L(t), 
cannot be determined without a precise knowledge of the microscopic dissipation mechanisms. 
However, since the dissipation of the magnetic energy involves the whole magnetosphere, it turns 
out that we may accurately reproduce the time variation of L(f) without a precise knowledge of 
the microscopic dissipative mechanisms. Indeed, in Sect. |5] we develop an effective description 
where our ignorance on the microscopic dissipative processes is encoded in few macroscopic pa- 
rameters, which allows us to determine the hght curves. In the remaining of the present Section 
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we investigate the spectral properties of the luminosity. To this end, we need to consider the syn- 
chrotron radiation spectral distribution. Since radiation from electrons is far more important than 
from protons, in the following we shall focus on electrons. It is well known that the synchrotron 
radiation will be mainly at the frequencies (IWallace 1977l|Grnzburg 1979^ : 

, eB 

CJm 7^ — , {11) 

trie 

where y is the electron Lorentz factor. Using Eq. M'il we get: 

co(r) ^ — {-] , R < r < lOR . (78) 



It is useful to numerically estimate the involved frequencies. To this end, we consider the giant 
flare of 1998, August 27 from SGR 1900+14: 

Bs ^ lA 10"* Gauss , — ^ 10"^ . (79) 

So that, from Eq. J78t it follows: 

(R \^ 
-j , R < r < \QR , (80) 

or 

wi ^ j^'?,.^! KeV < 0) < W2 - y^%.61 MeV . (81) 

The power injected into the magnetosphere is supplied by the azimuthal electric field during the 
initial hard spike. So that to estimate the total power we need to evaluate the power supplied by 
the azimuthal electric field. Let us consider the infinitesimal volume dV - sin OdrdOdip; the 
power supplied by the induced electric field in dV: 

6Bs ^ , 

dWE^ - riee ^ R^ sin^ 0drd0d<p , (82) 

spike 

where n^, is the electron number density. Since the magnetosphere is axially symmetric it follows 
that He cannot depend on (p. Moreover, within our theoretical uncertainties we may neglect the 
dependence on 6. So that, integrating over and tp we get: 

dWE^ ^ In^ n^e R^ dr . (83) 

spike 

In order to determine the spectral distribution of the supplied power, we note that to a good 
approximation all the synchrotron radiation is emitted at w,,,, Eq. J77> . So that, we may use 
Eq. (|80j to get: 

eB,\'" 1 



dr ^ - r \ — \ —du). (84) 

W3 



Inserting Eq. ( IMl into Eq. ( 1^ we obtain the spectral power: 

2n^ 6Bs „4 2 leBsV'^ 1 



F(lo) du ^ —r-ni,e v^R — doj , (85) 

i Ot spike \ ftle I (jji 



while the total luminosity is given by: 



J null 



L = F{(jj)du . (86) 
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Note that L is the total luminosity injected into the magnetosphere during the initial hard spike. 
So that, since the spike lasts 6 tspiUe, we have: 



Eburst - 5 hpike L , (87) 

where E hurst is the total burst energy. In the case of the 1998 August 27 giant burst from SGR 
1900+14 the burst energy is given by Eq. (I46> . Thus, using Eqs. (I87> and (I72> we have: 

L ^ 10« ^JIl , (88) 
sec 

which, indeed, is in agreement with observations. It is worthwhile to estimate the electron number 
density needed to dissipate the magnetic energy injected in the magnetosphere. To this end, we 
assume an uniform number density. Thus, using Eqs. (I86> . ( I85> and ( 18 H we get: 

L ^ \%T? n,e R'^ , (89) 

spike 

where we used ^ 1 . Specializing to the August 27 giant burst we find: 

Ue ^ 2.0 10'"* cm-^ , (90) 

indeed quite a reasonable value. Soon after the initial spike, the induced azimuthal electric field 
vanishes and the luminosity decreases due to dissipative processes in the magnetosphere. As 
thoroughly discussed in Sect. |5j it is remarkable that the fading of the luminosity can be accu- 
rately reproduced without a precise knowledge of the microscopic dissipative mechanisms. So 
that combining the time evolution of the luminosity L(t), discussed in Sect. [S] with the spectral 
decomposition we may obtain the time evolution of the spectral components. In particular, firstly 
we show that starting from Eq. ( I85> the spectral luminosities can be accounted for by two black- 
bodies and a power law. After that, we discuss the time evolution of the three different spectral 
components. 

The spectral decomposition Eq. ( I85t seems to indicate that the synchrotron radiation follows a 
power law distribution. However, one should take care of reprocessing effects which redistribute 
the spectral distribution. To see this, we note that photons with energy a> >2 nie quickly will pro- 
duce copiously almost relativistic pairs. Now, following Duncan & Thompson (1995), even if 
the particles are injected steadily in a time dtspike, it is easy to argue that the energy of relativistic 
particles is rapidly converted due to comptonization to thermal photon-pair plasma. Since the 
pair production is quite close to the stellar surface, we may adopt the rather crude approxima- 
tion of an uniform magnetic filed B ^ Bs throughout the volume Vpiasma - 12 nR^ (Duncan 
& Thompson 1995). Since typical magnetic fields in magnetars are well above Bqed, electrons 
and positrons sit in the lowest Landau levels. In this approximation we deal with an almost one 
dimensional pair plasma whose energy density is ( |Duncan & Thompson 1995) 1: 

_i 

2 2 I '^plasma \~ . fTlg 



Ue ^ niein,* +n,-) =^ j eBs '■ — exp(-- ) , (91) 

Cln) 2 \ / plasma 
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for Tpiasma < T plasma being the plasma temperature. Thus, the total energy density of the 
thermal photon-pair plasma is; 

4 

U = U, + Uy ^ Ue + — Tpi^^,^ . (92) 

The plasma temperature is determined by equaling the thermal energy Eq. (I92> with the fraction 
of burst energy released in the spectral region a> >2 nie. It is easy to find; 

Epairs - 0.147 £i„„, , (93) 

where for the numerical estimate we approximated wi - lO/TeVandwa - 10 MeV, correspond- 
ing to mildly relativistic electrons in the magnetosphere. So that we have; 



zr-i K \-^\ exp(-— ^) + - . ^ . (94) 

yLTl)- \ "^e I ^ plasma 'plasma 

In the case of August 27 giant burst from SGR 1900+14 we find: 

VI exp(-i) + 0.311 - 1.32 10"^ , x = , (95) 

X Me 

whose solution gives Tpiasma - 135 KeV. However, this is not the end of the whole story. Indeed, 
our thermal photon-pair plasma at temperature Tpiasma will be reprocessed by thermal electrons 
on the surface which are at temperature of the thermal quiescent emission Tq < \ KeV. So that, 
photons at temperature Tpiasma » Tq are rapidly cooled by Thompson scattering off electrons in 
the stellar atmosphere, which extends over several hundreds /erm/i beyond the edge of the star. 
The rate of change of the radiation energy density is given by (for instance, see lPeebles 1993> ; 

1 duy 4 CTj hq 



My dt Ifle 



(Tq - Tpiasma) , (96) 



where hq is the number density of electrons in the stellar atmosphere. The electron number 
density in the atmosphere of a p-star is of the same order as in strange stars, where hq ^ 10^^ cm'^ 
( |Alcock, Farhi, & Olinto 19'86t . So that, due to the very high electron density of electrons near 
the surface of the star, the thermal photon-pair plasma is efficiently cooled to a final temperature 
much smaller than Tpiasma- At the same time, the energy transferred to the stellar surface leads to 
an increase of the effective quiescent temperature. Therefore we are lead to conclude that during 
the burst activity the quiescent luminosity must increase. Let Ti be the final plasma temperature, 
then we see that the thermal photon-pair plasma contribution to the luminosity can be accounted 
for with an effective blackbody with temperature Ti and radius Ri of the order of the stellar 
radius. As a consequence the resulting blackbody luminosity is; 

Li = AjtRIctsbT^i , Ri < R . (97) 

In general, the estimate of the effective blackbody temperature Ti is quite difficult. However, 
according to Eq. ( I93> we known that Li must account for about 0.147 of the total luminosity. So 
that we have; 

Li(f) ^ 0. 141 Lit) . (98) 
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This last equation allow us to determine the blackbody temperature. For instance, soon after the 
hard spike we have L(0) - |^ - 10^^ ^ for the giant burst from SGR 1900+14. Thus, using 
R\-R, from Eq. we get; 

ri(0) ^ 61 KeV , (99) 

with surface luminosity Li(0) ^ lO'*^ 

Let us consider the remaining spectral power with w < Inie. We recall that the spectral power 
Eq. ( I85t originates from the power supplied by the induced electric field Eq. ( I83t . It is evident 
from Eq. ( I83> that, as long as ^ 1, the power supplied by the electric field does not depend 
on the mass of accelerated charges. Since the plasma in the magnetosphere is neutral, it follows 
that protons acquire the same energy as electrons. On the other hand, since the protons syn- 
chrotron frequencies are reduced by a factor the protons will emit synchrotron radiation near 
Wi. As a consequence, photons with frequencies near suffer resonant synchrotron scattering, 
which considerably redistribute the available energy over active modes. On the other hand, for 
w » the spectral power will follows the power law Eq. (I85> . Thus, we may write: 

Fioj) , 5wi<w<2m, , (100) 

where we have somewhat arbitrarily assumed the low energy cutoff ~ 5 <jJ\. On the other hand, 
for CO < 5 to\ the redistribution of the energy by resonant synchrotron scattering over electron 
and proton modes lead to an effective description of the relevant luminosity as thermal blackbody 
with eff'ective temperature and radius To and R2, respectively. Obviously, the blackbody radius 
R2 is fixed by the geometrical constrain that the radiation is emitted in the magnetosphere at 
distances r < 10 /J. So that we have: 

Ri < lOR . (101) 

The effective blackbody temperature T2 can be estimate by observing that the integral of the 
spectral power up to 5 u\ account for about the 60 % of the total luminosity. Thus, we have: 

L2{t) ^ 0.60 L(f) , (102) 

where 

L2 = AnRlo-sBT^2 ' Ri < lOR . (103) 

Equations ( I102> and ( I103> can be used to to determine the effective blackbody temperature. If we 
consider again the giant burst from SGR 1900+14, soon after the hard spike, assuming R2 - IQR, 
we readily obtain: 

7^2(0) ^ 21 KeV . (104) 

To summarize, we have found that the spectral luminosities can be accounted for by two black- 
bodies and a power law. In particular for the blackbody components we have: 

Ri < R , R2< lOR ; T2 < , Ri < R2 

(105) 

Li ^ 0.15 L , L2 - 0.60 L . 
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Interestingly enough, Eq. ( I105> displays an anticorrelation between blackbody radii and temper- 
atures, in fair agreement with observations. Moreover, the remaining 25% of the total luminosity 
is accounted for by a power law leading to the high energy tail of the spectral flux: 

dN 

, a ^ 2.33 , (106) 

dE 

extending up to £ ^ Im^ ^ 1 MeV. Indeed, the high energy power law tail is clearly displayed 
in the giant flare from SGR 1900+14 (see Fig. 3 in lFeroci et al. 1999t . and in the recent gigantic 
flare from SGR 1806-20 (see Fig. 4 in |Hurley et al. 2005| l. 

It is customary to fit the spectra with the sum of a power law and an optically thin thermal 
bremsstrahlung. It should be stressed that the optically thin thermal bremsstrahlung model is 
purely phenomenological. In view of this, a direct comparison of our proposal with data is 
problematic. Fortunately, Feroci et al. (2001) tested several spectral functions to the observed 
spectrum in the afterglow of the giant outburst from SGR 1900+14. In particular they found 
that, in the time interval 68 iec < t < 195 ^ec, the minimum spectral model were com- 
posed by two blackbody laws plus a power law. By fitting the time averaged spectra they re- 
ported dFeroci et al. 200 U : 

T2 ^ 9.3 KeV , 20.3 KeV , a 2.8 . (107) 

Moreover, it turns out that the power law accounts for approximately 10% of the total energy 
above 25 KeV, while the low temperature blackbody component accounts for about 85% of the 
total energy above 25 KeV. In view of our neglecting the contribution to energy from protons, we 
see that our proposal is in accordance with the observed energy balance. Unfortunately, in Feroci 
et al. (2001) the blackbody radii are not reported. To compare our estimate of the blackbody 
temperatures with the fitted values in Eq. J107> . we note that our values reported in Eqs. i99\ and 
ilQ4\ correspond to the blackbody temperatures soon after the initial hard spike. Thus, we need 
to determine how the blackbody temperatures evolve with time. To this end, we already argued 
that soon after the initial spike the luminosity decreases due to dissipative processes in the mag- 
netosphere. In Sect. |5lwe show that the fading of the luminosity can be accurately reproduced 
without a precise knowledge of the microscopic dissipative mechanisms. In particular, the rele- 
vant hght curve is given by Eqs. ( I123> and ( I127t . At f = we have seen that the total luminosity 
is well described by three different spectral components. During the fading of the luminosity, it 
could happens that microscopic dissipative processes modify the different spectral components. 
However, it is easy to argue that this does not happens. The crucial point is that the three spectral 
components originate from emission by a macroscopic part of the magnetosphere; moreover the 
time needed to modify a large volume of magnetosphere by microscopic processes is much larger 
than the dissipation time t^is ~ 10^ sec. Then we conclude that, even during the fading of the 
luminosity, the decomposition of the luminosity into three different spectral components retain 
its validity. Now, using the results in Sect.ls] we find: 

. 3.67 10- , ^^^i^J^ . 1.67 10- . (108) 
L(0) E({j) 
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Combining Eq. ( fTMt with Eqs. iTU% and ( fTU^t we obtain: 

T2(t - 68 sec) ^ ll.SKeV , Tiit - 195 sec) ^ 9.1 KeV , 
Tiit ^ 68 sec) ^ 26.1 KeV , T^{t ^ 195 sec) ^ 21.9 /TeV, 



(109) 



in reasonable agreement with Eq. MQll . Finally, let us comment on the time evolution of the 
spectral exponent a in the power law Eq. ( I106> . From Eq. ( I85> it follows that high energy modes 
have less energy to dissipate. Accordingly, once a finite amount of energy is stored into the 
magnetosphere, the modes with higher energy become inactive before the lower energy modes. 
As a consequence, the effective spectral exponent will increase with time and the high energy tail 
of the emission spectrum becomes softer. This explains also why the fitted spectral exponent a 
in Eq. (I107> is sUghtly higher than our estimate in Eq. ( I106> . 

4.4. HARDNESS RATIO 
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Fig. 4. Hardness-intensity plot of the time resolved hardness ratio, Eq. (II 10> . Data have been 
extract from Fig. 3 in Gotz et al. (2004). Full lines are our Eqs. il 15> and il 17> . 



Recently, it has been reported evidence for a hardness-intensity anti correlation within bursts 
from SGR 1806-20 (IGotz et al. 2004> . Indeed, Gotz et al. (2004) reported observations of the 
soft gamma ray repeaters SGR 1806-20 obtained in October 2003, during a period of bursting 
activity. They found that some bursts showed a significant spectral evolution. However, in the 
present Section we focus on the remarkable correlation between hardness ratio and count rate. 
Following Gotz et al. (2004) we define the hardness ratio as; 
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where H and S are the background subtracted counts in the ranges 40- 100 KeV and 15-40 KeV 
respectively. In Figure 0] we report the hardness ratio data extracted from Fig. 3 in Gotz et al. 
(2004). 

A few comments are in order. First, the hardness ratio becomes negative for large enough burst 
intensities. Moreover, there is a clear decrease of the hardness ratio with increasing burst in- 
tensities. Now we show that within our approach we may explain why the hardness ratio is 
negative and decreases with increasing burst intensities. To see this, we note that the hardness 
ratio Eq. ( II 10> is defined in terms of total luminosities in the relevant spectral intervals. Thus, to 
determine the total luminosity in the spectral interval ui - 1^2 we may use: 



•J ID 



L{a)i -W2) = F(u)dcD , (111) 



where F{u) is given by Eq. ( I85> . A straightforward integration gives: 

6Bs 

t spike 

where r\ and r2 are given by: 



Ucji - W2) - 2 7? He e - ^) , (112) 

Ots„ike \K Kl 



..,2.r^f^ A, . (113) 

me \r\,2i 



Assuming y ~ 1, we may rewrite Eq. ( II 13> as: 

wi,2 ^ lOMeV . (114) 

Using Eqs. (II 12> and (II 14> it is easy to determine the hardness ratio: 

HR - ^(^0 - - ^^^^ - ^ 1.66 - 2.44 ^ _ „ 

L(40- lOO/TeV) + L(15 -40/:ey) 1.66 + 2.44 

In Figure 0] we display our estimate of the hardness ratio Eq. dl 15> . We see that data are in 

quite good agreement with Eq. ( II 15> at least up to count rate ~ 5 lO-' counts j sec. For larger 

count rates data seem to lie below our value. We believe that, within our approach, there is a 

natural explanation for this effect. Indeed, for increasing count rates we expect that the hard tail 

w > 2 ^ 1 MeV of the spectrum will begin to contribute to the luminosity. According to 

the discussion in Sect. l4.3l these hard photons are reprocessed leading to an effective blackbody 

with temperature T\ . Now, for small and intermediate bursts the blackbody temperature T\ is 

considerably smaller than Eq. ( I99> . so that the effective blackbody contributes mainly to the soft 

tail of the spectrum. Obviously, the total luminosity of the effective blackbody is: 

L{\ - \0 MeV) =^ 2 jf- lie e R'^ {2.15 - \) ■ (116) 

spike 



Since this luminosity contributes to the soft part of the emission spectrum, Eq. (II 15> gets modi- 
fied as: 

1.66 - 3.59 

^^^L6^Ty59 ''''^ 
Equation (^^J is displayed in Fig. |3 for rates > 5 10^ counts I sec. Note that we did not take 
into account the proton contribution to the luminosity. Observing that protons contribute mainly 
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to luminosities at low energy o) < IQKeV, we see that adding the proton contributions leads to 
smaller hardness ratios bringing our estimates to a better agreement with data. In any case, we 
see that our theory allows to explain in a natural way the anti correlation between hardness ratio 
and intensity. 



5. LIGHT CURVES 

In our magnetar theory the observed burst activities are triggered by glitches which inject mag- 
netic energy into the magnetosphere where, as discussed in previous Sections, it is dissipated. As 
a consequence the observed luminosity is time depended. In this Section we set up an effective 
description which allows us to determine the light curves, i.e. the time dependence of the lumi- 
nosity. In general, the energy injected into the magnetosphere after the glitch decreases due to 
dissipative effects described in Sect. 14.31 leading to the luminosity L(t) = -^^^fp- Actually, the 
precise behavior of L(t) is determined once the dissipation mechanisms are known. However, 
since the dissipation of the magnetic energy involves the whole magnetosphere, we may ac- 
curately reproduce the time variation of L(t) without a precise knowledge of the microscopic 
dissipative mechanisms. Indeed, on general grounds we have that the dissipated energy is given 
by: 

L(t) ^ _dE(t)_ ^ ^(f-)^?, , 77 < 1 , (118) 
dt 

where 77 is the efficiency coefficient. Obviously the parameter K(t) does depend on the physical 
parameters of the magnetosphere. For an ideal system, where the initial injected energy is huge, 
the linear regime, where 77 = 1, is appropiate. Moreover, we expect that the dissipation time 
~ i is much smaller than the characteristic time needed to macroscopic modifications of the 
magnetosphere. Thus, we may safely assume K{t) ^ kq constant. So that we get: 

dE(t) 

m = - — ^ ^ K^E . (119) 
dt 

It is straightforward to solve Eq. ( II 19l i: 

E(t) = £0 exp(--) , L(f) = Lo exp(--) , ^ — , tq ^ - . (120) 

To To To /Co 

Note that the dissipation time tq = ^ encodes all the physical information on the microscopic 
dissipative phenomena. Since the injected energy is finite, the dissipation of energy degrades 
with the decreasing of the available energy. Thus, the relevant equation is Eq. ( II 18> with 77 < 1 . 
In this case, solving Eq. dl 18> we find: 

W) = Lo (l-— 1'" , (121) 



tdis 



where we have introduced the dissipation time: 



1 L/r" 

tdis = - 7^ . (122) 
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Then, we see that the time evolution of the luminosity is linear up to some time threak, after that 
we have a break from the linear regime rj - 1 to a non Unear regime with t] < 1 . If we indicate 
with tdis the total dissipation time, we get: 

t 

Lit) = Lo exp( ) , < f < tbreak , 

To 

(123) 



L(t) = Lithreak) I 1 - " " I , threak < t < tdis ■ 

Equation J123> is relevant to describe the light curve after a giant burst, where there is a huge 
amount of magnetic energy dissipated into the magnetosphere. It is interesting to compare our 
light curves, Eq. ( I123> . with the standard magnetar model. The decay of the luminosity in 
the standard magnetar model is due to the evaporation by a fireball formed after a giant burst 
and trapped onto the stellar surface ( |Duncan & Thompson I996| [Duncan & Thompson 200 1) . 
Indeed, Feroci et al. (2001) considered the light curves after the giant flare of 1998, August 27 
from SGR 1900+14, and the giant flare of 1979 March 5 from SGR 0526-66. Assuming fliat the 
luminosity varies as a power of the remaining fireball energy L ~ E", they found: 

L{t) = Lo , (124) 



where tevap is the time at which the fireball evaporates, and the index a accounts for the geom- 
etry and the temperature distribution of the trapped fireball. For a spherical fireball of uniform 
temperature a = |, so that the index a must satisfies the constrain: 

a <\ ■ (125) 

Note that our Eq. (I123> reduces to Eq. J124> if threak - and tj - a. However, we stress that our 
efficiency exponent must satisfy the milder constraint 77 < 1 . 

Feroci et al. (2001) performed a best fit of the light curve of the August 27 flare, background 
subtracted and binned to 5 sec, to Eq. (I124t and found: 

a = 0.756 ± 0.003 , t^^^p = 414 sec . (126) 

Indeed, from Fig. 2 in Feroci et al. (2001) one sees that the trapped fireball light curve account 
for the decay trend of the experimental light curve and matches the sudden final drop of the flux. 
However, it should stressed that the fitting parameter a in Eq. ( I126> does not satisfy the physical 
constraint Eq. jl25t . Even more, any deviations from spherical geometry or uniform temperature 
distribution lead to parameters a smaller than the upper bound |. Moreover, the trapped fireball 
light curve underestimate by about an order of magnitude the measured flux during the first stage 
of the outburst. We interpreted the different behavior of the flux during the initial phase of the 
outburst as a clear indication of the linear regime described by our Eq. ( I120> . As a matter of fact, 
we find that the measured light curve could be better described by Eq. (I123> with parameters (see 
Fig.E}: 

To = 8.80 sec , threak - 20 scc , T] - 0.756 , tdis - 414 sec . (127) 
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Fig. 5. Light curve after the giant flare of 1998, August 27 from SGR 1900+14. Dashed line is 
flie light curve in the standard magnetar model, Eq. J124> with parameters given in Eq. ( I126> . 
Full line is our light curve Eq. J123> with parameters in Eq. (I127> 

The same criticisms apply to the fit within the standard magnetar model of the light curve after 
the giant flare of 1979 March 5 from SGR 0526-66. The trapped fireball light curve fit in Feroci 
et al. (2001) gives: 

a = 0.71 ± 0.01 , tevap = 163 ± 5 sec . (128) 

Again the parameter a exceeds the bound Eq. jl25> . and the fit underestimates the flux during the 
first stage of the outburst (see Fig. 14 in lFeroci et al. 20()ll . Fitting our Eq. ( I123> to the measured 
flux reported in Fig. 14 in Feroci et al. (2001), we estimate: 

To = 15 sec , tbreak = 20 sec , Tj - 0.71 , tdis = 163 sec . (129) 

In Figs. |5l and |6l we compare our light curves Eqs. J123> . (I127> and ( I129> with the best fits 
performed in Feroci et al. (2001). Obviously, both Ught curves agree for t > tbreak, while in the 
linear regime t < tbreak, where the trapped firebafl light curves underestimate the flux, our light 
curves follow the exponential decay and seem to be in closer agreement with observational data. 

Several observations indicate that after a giant burst there are smaller and recurrent bursts. 
According to our theory these small and recurrent bursts are the effect of several small glitches 
foUowing the giant glitch. We may think about these small bursts like the seismic activity fol- 
lowing a giant earthquake. These seismic glitches are characterized by Ught curves very diff'erent 
from the giant burst light curves. In the standard magnetar model these light curves are accounted 
for with an approximate f "-^ decay ( |Lyubarski, Eichler, & Thompson 2002) . Within our theory 
there is a natural way to describe the seismic burst activity. Indeed, during these seismic bursts, 
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Fig. 6. Light curve after the giant flare of 1979 March 5 from SGR 0526-66. Dashed line is the 
light curve in the standard magnetar model, Eq. (I124t with parameters given in Eq. ( I128> . Full 
line is our light curve Eq. (I123> with parameters in Eq. il2% 



that we shall call settling bursts, there is an almost continuous injection of energy into the mag- 
netosphere which tends to sustain an almost constant luminosity. This corresponds to an effective 
K in Eq. Jl 18t which decreases smoothly with time. The simplest choice is: 



Inserting into Eq. ( II 18> and integrating, we get: 
E(t) = 

So that the luminosity is: 



1 + Kit 



eI" -(1-7?)- ln(l+^if) 



Lit) 



(l+Klt)" 

After defining the dissipation time: 



1 -(1-77) 



KlE. 



ln(l + Kit) 



ln(l +Kitdi,) = 



'^1 ^0 



/Co 1 



: rewrite Eq. ([22) ■ 



{\+Klt)1 



1 - 



ln(l + Kit) 



ln(l +Kitdis) 



(130) 



(131) 



(132) 



(133) 



(134) 



Note that the light curve Eq. ( I134t depends on two characteristic time constants ^ and f^/j. We 
see that k\ f^/j, which is roughly the number of small bursts occurred in the given event, gives an 
estimation of the seismic burst intensity. Moreover, since during the seismic bursts the injected 
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energy is much smaller than in giant bursts, we expect that fitting Eq. jl34t to the observed light 
curves will result in values of rj smaller than the typical values in giant bursts. In the following 
Sections we show that, indeed, our light curves Eq. J134l i are in good agreement with several 
observations. 

5.1. AXP 1E 2259+586 

On 2002, June 18 SGR-like bursts was recorded from AXP IE 2259+586. Coincident with the 
burst activity were gross changes in the pulsed flux, persistent flux, energy spectrum, pulse pro- 
file and spin down of the source JWoods et al. 2004b . As discussed in previous Sections, these 
features are naturally accounted for within our theory. However, we believe that the most remark- 
able and compelUng evidence for our proposal comes from the observed coincidence of the burst 
activity with a large glitch. Moreover, the time evolution of the unabsorbed flux from AXP IE 
2259+586 following the 2002 June outburst reported in Woods et al. (2004) can be explained nat- 
urally within our theory. We could consider the June 18 SGR-like bursts from AXP IE 2259+586 
the Rosetta Stone for our magnetar theory. 

The temporal decay of the flux during the burst activity displays a rapid initial decay which lasted 
about 1 days, foUowed by a more mild decay during the year following the onset of the burst ac- 
tivity. Indeed, Woods et al. (2004) splitted the data into two segments, and fit each independently 
to a power law: 

F(t) ~ f«> , ai = - 4.8 ±0.5 , f < 1 days , 

(135) 

F(t) ~ f',a2 = - 0.22 ± 0.01 , f > 1 days . 

It is evident from Eq. ( I135> that the standard magnetar model is completely unable to reproduce 
the phenomelogical power law fit. On the other hand, even the phenomenological parametrization 
cannot account for the time evolution of the flux. Indeed, if we assume the power law Eq. (I135> 
for the decay of the flux, then we cannot explain why and how the source returns in its quiescent 
state with quiescent flux JWoods et al. 20041 : 

Fn ^ 1.53 10"" _f^£f_ (j3g-) 

cm-^ sec 

Note that adding the quiescent flux to the power law decay does not resolve the problem, for in 
that case the fits worst considerably. Our interpretation of the puzzling light curve displayed in 
FigQis that AXP IE 2259+586 has undergone a giant burst at the glitch epoch, and soon after 
the pulsar has entered into a intense seismic burst activity. Accordingly, the flux can be written 
as: 

Fit) = FcBit) + FsBit) + Fq , (137) 

where Fq is the quiescent flux, Eq. ( I136> , Fcgit) is the giant burst contribution to the flux given 
by Eq. jl23> . and Fs sit) is the seismic burst contribution given by Eq. jl34> . Since during the first 
stage of the outburst there are no available data, we may parameterize the giant burst contribution 
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Fig. 7. The time evolution of the unabsorbed flux from AXP IE 2259+586 foUowing the 2002 
June outburst. Data have been extracted from Fig. 13 in Woods et al. (2004). Dashed lines are the 
phenomelogical power law fits, Eq. ( I135> . Full line is our light curve Eq. ( I137t . with parameters 
in Eq. iTM . 



as: 

FcB(t) = FcBiO) 1 1 
while Fs sit) is given by: 

FsB(t) = 



'IGB 



Fsb(0) 

(1 +Klt)'T' 



1 - 



t 

tcB 



ln(l + Kit) 



< t < tcB 



Q < t < tsi 



(138) 



(139) 



ln(l + KytsB) _ 

where tcB and tss are the dissipation time for giant and seismic bursts respectively. In Fig.Qwe 
display our light curve Eq. J137l i with the following parameters: 



FcBiO) ^ 1.5 10"** _fl£i_ , T^^g ^ 0.828 , tcB ^ 0.91 days 



Fsb(O) - 5.0 10 



cm^ sec 
-11 ergs 



(140) 



VSB 



0.45, tsB - 10 days , ki ^ 0.20 days 



cm-^ sec 

A few comments are in order. First, the agreement with data is rather good. Second, our efficiency 
exponent ricB is consistent with the values found in the giant bursts from SGR 1900+14 and SGR 
0526-66. On the other hand, quite consistently, we have tjsb < J]gb- Finally, we stress that from 
our interpretation of the light curve it follows that the onset of the intense seismic burst activity 
{K\tsB ~ 200) did not allow a reliable estimation of ^, which we predicted to be of order 10"^. 
Interestingly enough, following the 2002 June outburst it was detected an infrared flux change 
correlated with the X-ray flux variability ( ITam et al. 2004> . Since the observations began three 
days after the 2002 June outburst, according to our theory the infrared flux is parameterized as: 



FiVt) 



(1 + Klt)''^" 



1 



ln(l + Kit) 



In(l + Kitss) 



'ISB 



+ F'; 



< t < tsi 



(141) 
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Fig. 8. The time evolution of the unabsorbed IR flux from AXP IE 2259+586 following the 
2002 June outburst. Data have been extracted from Fig. 1 in Tam et al. (2004). Dashed line is the 
phenomelogical power law fit f-0-2i±o.oi p^jj Yine is our light curve Eq. J141> . 

with the same parameters as in Eq. jl39> . Indeed, assuming F'J^^iQ) ^ 9.5 ^^^^^ cnr-^Tec-' ^^'^ 
Fq ^3.3 10~'^ cwrVc-i ' found that our light curve Eq. MAW is in remarkable good agreement 
with data (see Fig. The strong correlation between infrared and X-ray flux decays observed 
after the 2002 June outburst from AXP IE 2259+586 strongly suggests a physical link between 
the origin of both type of radiation. 

5.2. SGR 1900+14 

Soon after the 1998, August 27 giant burst, the soft gamma repeater SGR 1900+14 entered a 
remarkable phase of activity. On August 29 an unusual burst from SGR 1900+14 was detected 
JIbrahim et al. 200'll which lasted for a long time ~ lO^iec. As discussed in Ibrahim et al. (2001), 
on observational grounds it can be ruled out extended afterglow tails following ordinary bursts. 
In Figure |9] we display the flux decay after the August 29 burst. Data has been extracted from 
Ibrahim et al. (2001). In Ibrahim et al. (2001) the temporal behavior of the flux decay has been 
parameterized as a power law (dashed line in Fig.|5): 

Fit) = (89.16 + 1.34) 10-'°^^ r<»*'™25). (142) 

aw^ sec 

As already stressed, the phenomenological power law decay cannot explain the return of the 
source in its quiescent state where the flux is ( [Hurley et al. 1999| l: 

Fn = 0.96 ±0.07 10'" ^'y^ . (143) 

cm^ sec 
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Fig. 9. Flux evolution of the August 29 burst from SGR 1900+14. Data has been extracted from 
Fig. 4, panel (d), in Ibrahim et al. (2001). Dashed line is the phenomenological power law fit 
Eq. J142> : full line is our Ught curve Eqs. (I144> and ( I145> . 
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Fig. 10. The time evolution of the unabsorbed flux from SGR 1900+14 following the 1998 
August outburst. Data has been extracted from Fig. 2 in Woods et al. (2001). Dashed line is the 
power law best fit F(t) ~ r(°-^'3±° °25) line is our Hght curve Eqs. J144l i and ( I146I I. 



On the other hand, we may easily account for the observed flux decay by our light curve: 



F(t) 



F(0) 

(1 + Kit)" 



ln(l + Kit) 



ln(l + Kitdis) 



+ Fr 



(144) 
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where Fq is fixed by Eq. ( I143> . Indeed, in Fig.|5]we compare our light curve Eq. J144l i with 
observational data. The agreement is quite satisfying if we take: 

F(0) ^ 1.05 10"'^ _^!;gf_ ^ ri ^ Q.5, tdi, ^ 1.2 10^ sec , ki ^ 0.50 sec'' . (145) 
cm^ sec 

Woods et al. (2001) have analyzed a large set of X-ray observations of SGR 1900+ 14 in order 
to construct a more complete flux history. They found that the flux level was more than an order 
of magnitude brighter than the level during quiescence. This transient flux enhancement lasts 
about 40 days after the giant flare. Unlike Woods et al. (2001), that argued that this enhancement 
was an artifact of the August 27 flare, we believe that the flux history can be adequately described 
as seismic burst activity of the source. In Fig.[TO|we report the flux light curve extracted from 
Fig.2 in Woods et al. (2001) together with their power law best fit. Again we find the the flux 
history is accounted for quite well by our light curve Eq. ( I144t with the following parameters: 

F(0) ^ 4.8 10"^ _f^£f_ ^ ^ ^ Q 55 ^ ^^.^ ^ 2OO days , ki ^ 2 10^ days'^ . (146) 
cm-^ sec 

The agreement between our light curve Eqs. J144t and ( I146t with the power law best fit is strik- 
ing. Moreover we see that our curve deviates from the power law fit for f > 60 days tending to Fq 
at f = tdis- Woods et al. (2001) noted that extrapolating the fit to the August 27 X-ray light curve 
back toward the flare itself, one finds that the expected flux level lies below the ASM flux mea- 
surements (squares in Fig.llO>. Moreover, these authors observed that the discrepancy reduces 
somewhat when one pushes forward the reference epoch to about 14 minutes after the onset of 
the flare. We believe that the discrepancy is due to a true physical effect, namely the observed 
discrepancy from extrapolated light curve and ASM measurements is a clear indication that the 
surface luminosity increases after the burst activity. In particular, soon after the August 27 giant 
flare we have seen in Sect. 14.31 that the surface temperature increases up to ~ 61 KeV and the 
surface luminosity reaches ~ lO'*'* Almost all the deposited energy is dissipated within the 
dissipation time of the giant flare ~ 400 sec. Nevertheless, it is natural to expect a more gradual 
afterglow where a small fraction of the energy deposited onto the star surface is gradually dissi- 
pated. As a matter of fact, we find that the observed level of luminosity Lx ~ 10^^ ^ (assuming 
a distance d - 10 kpc) at about 0.01 days since the August 27 giant flare, is consistent with the 
gradual afterglow scenario. 

On 2001 April 18 the soft gamma ray repeater SGR 1900-H14 emitted an intermediate burst. 
The light curve of this event did not show any initial hard spike and was clearly spin-modulated. 
Moreover, the energetics appeared to be intermediate in the 40 - 100 KeV range, with a total 
emitted energy of about 1.9 10^^ ergs dPeroci et al. 20031 . In Fig.E]we report the temporal be- 
havior of the X-ray (2 - 10 KeV) flux from SGR 1900+14 in the aftermath of flie 2001 April 18 
flare. Data has been extracted from Fig. 2 in Feroci et al. (2003). Feroci et al. (2003) attempted a 
simple power law function to the flux data: 



Fit) ~ r" + K 



(147) 
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Fig. 11. Temporal behavior of the X-ray flux from SGR 1900+14 in the aftermath of the 2001 
April 18 flare. Data have been extracted from Fig. 2 in Feroci et al. (2003). Dashed line is the 
power law best fit Eqs. ( I147> and (I148> . Full line is our light curve Eqs. (I144> and ( I149> . 

where the constant K should take care of the quiescent luminosity. Indeed, Feroci et al. (2003), 
fitting Eq. (I147> to the data, found: 



a = 0.89(6) , K = 0.78(5) 10" 



(148) 



cm^ sec 

As it is evident from Fig. [H] the power law globally fits the data quite nicely. However, the 
reduced turns out to be in excess to 3, mainly due to the bump occurring in the light curve 
at f ~ 10^ sec JFeroci et al. 2003t . Indeed, after excluding the bump they get a good fit with 
Idof ^ 1 without an appreciable variation of the fit parameters ( IFeroci et al. 2003l l. However, 
there is still a problem with the phenomelogical power law decay of the flux. As a matter of fact, 
Eq. ( I148> shows that the power law fit underestimates the quiescent luminosity. In our opinion 
this confirms that the phenomenological power law decay of the flux is not adeguate to describe 
the time variation of the flux. On the other hand, we find that our light curve Eq. ( I144> . with 
quiescent luminosity fixed to the observed value Eq. ( I143> . furnishes a rather good description 
of the flux decay once the parameters are given by: 



F(Q) ^ 2.6 10"^ _£^£i_ ^ ^ ^ Q gg ^ j^.^ ^ 3 jq6 ^ ^ q_25 sec'^ 



(149) 



cm'- sec 

Within our interpretation, Eq. MA91 shows that the flux decay in the aftermath of the April 18 
flare is characterized by a very large seismic burst activity (/ciff/,, ~ 10*), which lasts for about 
10* sec. So that the bump in the flux at f ~ 10^ sec is naturally explained as fluctuations in the 
intensity of the seismic bursts. 

Lenters et al. (2003) reported the spectral evolution and temporal decay of the X-ray tail of a 
burst from SGR 1900+14 recorded on 2001 April 28, 10 days after the intense April 18 event. In 
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time ( sec ) 

Fig. 12. The temporal decay of the flux from the 2001 April 28 burst from SGR 1900+14 in 
the energy band 2-20 KeV. The data has been extracted from Fig. 5 in Lenters et al. (2003). 
Dashed line is the phenomenological best-fit power law times exponential function adopted in 
Lenters et al. (2003) to describe data, Eqs. ( I150> and ( I151> . Full continuous line is our light curve 
Eqs. ([HU and ( fT52b . 

Fig-Elwe display the temporal decay of the flux from the 2001 April 28 burst in the energy band 
2-20 KeV. The data has been extracted from Fig. 5 in Lenters et al. These authors attempted 
several functional forms to fit the decay of the flux. They reported that the decay was equally 
wefl fitted by either a power law times exponential or a broken power law. We stress that both 
fits are phenomenological parametrization of the observational data, and that both fits are unable 
to recover the quiescent flux. For definitiveness, we shall compare our light curve with the power 
law times exponential fit: 

F(t) ~ r" exp(--) , (150) 

T 

The best fit to the temporal decay of the flux from the 2001 April 28 burst in the energy band 
2-20 KeV gives JLenters al. 2003> : 

a =0.68 + 0.04 , r = 5 ± 1 10^ sec . (151) 

In Figure El we compare the phenomenological best fit Eqs. ( I150> and ( I151> with our light 
curve Eq. J144> . where the quiescent luminosity is fixed to the observed value Eq. J143t . and the 
parameters are given by: 

F(0) ^ 1.4 10"'' _f^^ , 77 ^ 0.5 , tdis ^ 5.5 10^ sec , Ky ^ 0.06 sec'^ . (152) 
cm^ sec 

Again, we see that our light curve gives a quite satisfying description of the flux decay. 

Let us consider, finally, the light curve for the intermediate burst from SGR 1900+ 14 occurred on 
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Time (seconds) 



Fig. 13. Time history of the 2001 July 2 burst from SGR 1900+ 14 in the energy band 7-100 KeV 
as observed by FREGATE. The data has been extracted from Fig. 1 in Olive et al. (2004) by 
binning the light curve histogram. Full line is our light curve Eqs. ( I144> and (I153> . 

2001 July 2 JOlive et al. 2004> . In Figure [01 we display the time decay of the flux after the July 
2 burst. The data have been extracted from Fig. 1 in Olive et al. (2004) by binning the light curve 
histogram. The displayed errors are our estimate, so that the data should be considered as purely 
indicative of the decay of the flux. We find that Fig. 1 in Olive et al. (2004) is very suggestive, for 
it seems to indicate that the burst results from several small bursts, i.e. according to our theory 
the burst is a seismic burst. As a consequence we try the fit with our light curve Eq. J144> . In this 
case the quiescent flux has been fixed to Fq ^ 0, for the observational data has been taken in the 
energy range 7 - 100 KeV where the quiescent flux is very small. Attempting the fit to the data 
we find: 

A counts 1 

F{0) ^ 4.07 10^ , 77 ^ 0.36 , t^s ^ 40 sec , ki ^ 5.0 sec'^ . (153) 

sec 

The resulting light curve is displayed in Fig.^] The peculiarity of this burst resides in the fact 
that the burst activity terminates suddenly at f ^ 4 sec well before the natural end at fj, 5 ^ 40 sec. 

5.3. SGR 1627-41 

SGR 1627-41 was discovered with the Burst And Transient Source Experiment (BATSE) 
on the Compton Gamma-Ray Observatory (CGRO) in June 1998 ( IKouvehotou et al. 19981 
IWoods et al. 1999at when it emitted over 100 bursts within an interval of 6 weeks. Kouveliotou 
et al. (2003) presented the results of the monitoring of the flux decay of the X-ray counterpart of 
SGR 1627-41 spanning an interval of roughly five years. Moreover, these authors attempted to 
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Fig. 14. Time decay of the flux from SGR 1627-41. The data has been taken from Table 1 in 
Kouveliotou et al. (2003). Full line is our best fitted Hght curve Eqs. (I144> . ( I154> and ( I155> . 

understand within the standard magnetar model the three year monotonic decline of SGR 1627- 
41 as cooling after a single deep crustal heating event coinciding with the burst activity in 1998. 
They assumed an initial energy injection to the crust of the order of 10^* ergs. However, it must 
be pointed out that this assumption is highly unrealistic, for the estimate of the total energy re- 
leased in bursts during the activation of SGR 1627-41 ranges between 4 10"^^ - 2 lO^-' ergs. 
In addition, since gamma rays was not detected, they assumed that the conversion efficiency of 
the total energy released during the activation into soft gamma rays were considerable less than 
100 %. They also assumed that the core temperature is low, i.e. the core cools via the direct 
URCA process. Notwithstanding these rather ad hoc assumptions, Kouveliotou et al. (2003) was 
unable to explain the March 2003 data point, which clearly showed that the flux did not decay 
further (see Fig.ll4>. In other words, the levelling of the flux during the third year followed by 
its sharp decline is a feature that is challenging the standard magnetar model based on neutron 
stars, and that begs for an explanation within that model. On the other hand, we now show that 
the peculiar SGR 1627-41 light curve find a natural interpretation within our theory. In Fig. 1141 
we display the time decay of the flux. The data has been taken from Table 1 in Kouveliotou et al. 
(2003). In this case we are able to best fit our light curve Eq. (I144> to available data. Since the 
number of observations is rather low, to get a sensible fit we have fixed the dissipation time to 
1200 days and the quiescent luminosity to the levelling value at f > 1200 days: 



Fq ^ 2.7 10" 

cm^ sec 

The best fit of our light curve to data gives: 

^-11 ^''8^ 



tdis 



1200 days 



FiO) = 0.83(11) 10" 



cm-^ sec 



0.25(8) , Ki = 0.04(1) days- 



USA) 



(155) 
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with a reduced ~ 1. From Fig.^^ where we compare our best fitted light curve with data, 
we see that our theory allow a quite satisfying description of the three year monotonic decline of 
the flux. 



5.4. SGR 1806-20 

SGR 1806-20 entered an active phase in 2003, culminating in a gigantic flare on 2004 December 
27, with energy greatly exceeding that of all previous events. Figure 3 in Hurley et al. (2005) 
displays the time history of the flux averaged over the rotation period of the pulsar soon af- 
ter the giant flare. These authors fitted the light curve within the standard magnetar model 
based on the evaporation of a fireball formed after the giant burst and trapped onto the stel- 
lar surface, Eq. (I124t . The fit of the rotation smoothed curve to the fireball function gives 
( |Hurley et al.^005> : 

a = 0.606 ± 0.003 , f,™^ = 382 + 3 sec . (156) 
However, from Figure 3 in Hurley et al. (2005) it is evident that fireball light curve underestimate 

' \ ' 

9.04 - 
9 - 
8.96 - 

O " 

- 8.92 - 



a 




2000 4000 

Time (s) 

Fig. 15. Time history of the second component after the 2004 December 27 giant burst from 
SGR 1806-20. The data has been extracted from Fig. 5 in Mereghetti et al. (2005b) by binning 
the light curve histogram. The continuous curve is our light curve Eqs. J144> . ( I157> and ( I158t . 



the luminosity for f < 30 sec. Thus we see that the Ught curve can be better accounted for by our 
light curve Eq. ( I123> with tbreak - 30 sec, quite close to the values found for the giant bursts from 
SGR 1900+14 and SGR 0526-66. A more precise determination of the parameters of our light 
curve, however, must await for more precise data in the initial phase of the afterglow. Instead, 
in the present Section we discuss the light curve of a second, separate component after the giant 
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burst reported in Mereghetti et al. (2005b). Indeed, Mereghetti et al. (2005b) found evidence for 
a separate component in the light curve starting at t ~ 400 sec from the onset of the giant burst, 
forming a peak at f ~ 600 sec and ending at f ~ 3000 sec (see Fig. 5 in |Mereghetti et al 2005b| l. 
As already discussed, in our theory it is expected that there is an intense seismic burst activity 
following a giant burst. In Figure[2]we display the flux history starting from the giant flare. We 
show a few points of the second component extracted from Fig. 5 in Mereghetti et al. (2005b) by 
binning the light curve histogram. The displayed errors are our estimate, so that the data should 
be considered as purely indicative of the decay of the flux. We fit the data with our fight curve 
Eqs. lll44> assuming: 

4 count 

Fq ^ 8.796 10 , tdis = tend - tstart , tstart - 625 SeC , tend - 3000 SeC . 

sec 

(157) 

Tfie best fit of our light curve to data gives: 

A count , 

F(0) = 0.074 10^ , 1] = 0.18 , Ki = 0.10 sec'^ . (158) 

sec 

In Fig.^lwe compare our light curve with data and find that our theory allow a quite satisfying 
description of the time history of the flux. 

6. CONCLUSIONS 

Let us summarize the main results of the present paper We have discussed p-stars endowed 
with super strong dipolar magnetic field. We found a well defined criterion to distinguish ro- 
tation powered pulsars from magnetic powered pulsars (magnetars). We showed that glitches, 
that in our magnetars are triggered by magnetic dissipative efifects in the inner core, explain both 
the quiescent emission and bursts from soft gamma-ray repeaters and anomalous X-ray pulsars. 
In particular, we were able to account for the braking gUtch from SGR 1900H-14 and the nor- 
mal glitch from AXP IE 2259+586 following a giant burst. We accounted for the observed anti 
correlation between hardness ratio and intensity. Within our magnetar theory we were able to 
account quantitatively for light curves for both gamma-ray repeaters and anomalous X-ray pul- 
sars. In particular we explained the light curve after the June 18, 2002 giant burst from AXP IE 
2259+586. The abiUty of our p-star theory to reach a satisfying understanding of several observa- 
tional features of soft gamma-ray repeaters and anomalous X-ray pulsars supports our proposal 
for a revision of the standard paradigm of relativistic astrophysics. 

Let us conclude by briefly addressing the theoretical foundation of our theory. As a mat- 
ter of fact, our proposal for p-stars stems from recent numerical lattice results in QCD 
JCea & Cosmai 20031 . which suggested that the gauge system gets deconfined in strong enough 
chromomagnetic field. This leads us to consider the new class of compact quark stars made 
of almost massless deconfined up and down quarks immersed in a chromomagnetic field in (3- 
equilibrium. Our previous studies showed that these compact stars are more stable than neutron 
stars whatever the value of the chromomagnetic condensate. This remarkable result indicates that 
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the true ground state of QCD in strong enough gravitational field is not realized by hadronic mat- 
ter, but by p-matter In other words, the final collapse of an evolved massive star leads inevitably 
to the formation of a p-star. 
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